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The Spitzer-IRAC Point Source Catalog of the Vela-D Cloud 
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ABSTRACT 

This paper presents the observations of the Cloud D in the Vela Molecular 
Ridge, obtained with the IRAC camera onboard the Spitzer Space Telescope 
at the wavelengths A= 3.6, 4.5, 5.8, 8.0 /im. A photometric catalog of point 
sources, covering a field of approximately 1.2 square degrees, has been extracted 
and complemented with additional available observational data in the millimeter 
region. Previous observations of the same region, obtained with the Spitzer 
MIPS camera in the photometric bands at 24 fim and 70 fim, have also been 
reconsidered to allow an estimate of the spectral slope of the sources in a wider 
spectral range. A total of 170,299 point sources, detected at the 5 a sensitivity 
level in at least one of the IRAC bands, have been reported in the catalog. There 
were 8796 sources for which good quality photometry was obtained in all four 
IRAC bands. For this sample, a preliminary characterization of the young stellar 
population based on the determination of spectral slope is discussed; combining 
this with diagnostics in the color-magnitude and color-color diagrams, the relative 
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population of young stellar objects in the different evolutionary classes has been 
estimated and a total of 637 candidate YSOs have been selected. The main 
differences in their relative abundances have been highlighted and a brief account 
for their spatial distribution is given. The star formation rate has been also 
estimated and compared with the values derived for other star forming regions. 
Finally, an analysis of the spatial distribution of the sources by means of the 
two-point correlation function shows that the younger population, constituted 
by the Class I and flat-spectrum sources, is significantly more clustered than the 
Class II and III sources. 

Subject headings: infrared: stars — ISM: clouds — ISM: individual (Vela Molec- 
ular Ridge) — stars: formation — stars: pre- main- sequence — surveys 



Introduction 



The Vela Molecul ar Ridge (hereinafter V MR) is a cloud complex constituted by four 
giant molecular clouds (iMurphy fc May Ill99ll ). hosting large regions of a ctive star formatio n 
(hereinafter SF), in which low and intermediate mass stars are forming (ILiseau et al.lll992l ). 
Its position on the Galactic plane (272° > £ > 260°; 3° > 6 > —3°) makes this complex a 
unique laboratory for studying the SF phenomenon when the environmental conditions are 
those typical of the Galactic disk. This is an important point to consider in extrapolating 
global star formation parameters to the Galactic scale. Indeed, the most commonly studied 
SF regions are projected outside the Galactic plane, in part because of the observational 
need to minimize the source confusion. The drawback is that the global quantities derived 
for these regions could not be fully representative of the bulk of the SF in the Galaxy. 



Our study is focused on a p articular subregion of the Vela cloud complex, the so-called 
D cloud (IMurphy fc May Ill99ll . hereinafter VMR-D), that has already been the subject of 



previous papers by our group. Former s tudies have shown a rich phenomenology a ssoci- 



(Lorenzetti et al. 



(Massietal 



2000 



Yamaguchi et al. 



2002 



ated with this cloud, involving outflows (IWouterloot fc Brandl Il999t lElia et al.l 120071 ) , jets 



Giannini et al. 20051 ). and evidences of clustered and isolated SF 



20031) . After the CO emission survey of the whole VMR carried out by 
(Il999l ). two maps of ^ 1 square degree region of VMR-D were obtained 
with higher resolu tion by lElia et al.l (120071 ) (^^CO(l-O) and ^^C0(2-l) emission lines) and 



Massi et al.l (120071 ) (mm dust continuum), aimed at studying the cloud morphology and 
properties as well a s their relationships with the SF process active in this region (see also 
De Luca et al.ll2007l ). This molecular cloud has been recently observed with both MIPS and 
IRAC cameras onboard the Spitzer Space Telescope and the MIPS aperture photometry has 
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already been the subject of a previous work (iGiannini et al.l 120071 ) 



In this paper we present the analysis of the Spitzer-IRAC mosaics of the VMR-D cloud, 
covering a region which mostly overlaps the field observed at longer wavelengths by MIPS 
and also corresponding to the sky regions already mapped in the mm spectral range and 
discussed in the previous papers. These new IRAC observations allowed us to compile 
a four-band photometric catalog of point sources that, complemented by the information 
previously obtained at longer wavelengths, constitutes a unique data bank for studying the 
young stellar objects (YSOs) in this cloud. With this as our goal we have also revisited 
the photometric analysis of the MIPS 24 /xm and 70 fim mosaics to obtain more accurate 
PSF photometry. The inclusion of the longer wavelength photometry is motivated by our 
aim to determine, as accurately as possible, the spectral slopes of the sources, infer their 
evolutionary stages, and obtain a reliable statistical description of the YSO population. 

A particularly valuable characteristic of our IRAC observations is that the data were 
collected by splitting the total integration time into two observing epochs 4.5 months apart, 
thereby enabling a first exploration of source variability in VMR-D. Indeed we found a 
number of sources that clearly showed brightness variations; they constitu te a significant 



subsa mple deserving special attention, and are discussed in a separate paper (IGiannini et al. 



20091 ). 



This paper is organized as follows: after a brief description of the observations in § HI 
we give an account of the procedure followed for source extraction and photometry for both 
IRAC and MIPS data in § |3l In the same section we also discuss the source selection 
criteria as well as the completeness and reliability. The final VMR-D point source catalog is 
presented in §|11 In §|5]we focus our attention on a subsample of 8796 sources (out of a total 
of 170299) that have been detected in all the four IRAC bands and for which good quality 
photometry was obtained. In the following we shall refer to this subsample as the "working 
sample". By means of color-color and color-magnitude diagrams, we select the candidate 
YSOs obtaining a census of the YSO content of the VMR-D region. In the same section a 
brief account of the clustering properties is done. Our conclusions are summarized in § |6l 



Observations 



^The observations were carried out with the Infrared Array C amera (IRAC; iFazio et al. 

20041 ) onboard the Spitzer Space Telescope ( iWerner et al.l 120041 ). The data were obtained 



as part of the Spitzer/IRAC Cycle 3 GTO (PID 30335, PI G. Fazio). The Galactic region 
observed was dehmited in longitude by 264°29' >i> 263°00' and in latitude by 0°42' > 6 > 
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— 0°50', covering about 1.2 square degrees. This field corresponds to a significant portion 
of the VMR-D cloud that has been mapped in four spectral bands, corresponding to the 
four IRAC camera detectors observing at the wavelengths A = 3.6, 4.5, 5.8, and 8.0 fim, 
respectively. The effective integration time per pixel, texp = 20.8 s in each band, has been 
accumulated by splitting the observations in two epochs separated by 4.5 months, a choice 
mainly dictated by the desire to monitor the observed region for possible variability of the 
sources. The same splitting of the observation time into two epochs has also been beneficial 
for a better cleaning of artifacts from the final mosaics. 

The IRAC data were obtained by mapping the whole region, in each epoch, with 12x14 
tiles, each observed in two dither positions. The observation was performed in IRAC 12 s 
HDR mode, consisting of two consecutive integrations of 0.4 and 10.4 s exposure times. The 
HDR mode was chosen to limit the saturation effects around the bright regions, a problem 
which is clearly present in the long exposure images in the neighborhood of the brightest 
sources. The data were reduced using the IRACproc package (Schuster et al. 2006) to 
combine the individual tiles into a flux calibrated mosaic on a pixel grid of 0.8627" /pix. 
IRACproc is a wrapper for the Spitzer Science Center mosaic software MOPEX and provides 
enhanced outlier (cosmic rays) rejection. 

Thus for each spectral band we obtained a total of four mosaics, one long and one short 
exposure for each of the two epochs. The mosaics of each epoch were subsequently coadded to 
obtain two flnal mosaics in each IRAC band, one with the long [t^xp = 20.8 s) and the other 
with the short {t^xp = 0.8 s) total exposure time. However, due to the different orientation 
constraints in the two observing epochs, and to the fact that the four IRAC cameras do 
not see exactly the same flelds simultaneously, the images obtained in the two epochs do 
not fully overlap. As a consequence, in the coadded mosaic the region corresponding to 
the full exposure time is limited by the spatial overlap of the mosaics obtained in the two 
observing epochs. The flnal mosaics therefore include a central region, observed for the 
total integration time texp = 20.8 s, as well as adjacent surrounding regions corresponding to 
the non-overlapping flelds that have been imaged only in one of the two epochs for a total 
exposure time texp = 10.4 s. The geometry of the mapping and the different parts covered 
are shown in Figure [TJ In the following we shall consider the whole region observed in the 
four bands, delimited in Figured] by the yellow contour, to produce the corresponding point 
source catalog, recalling that the area that has been imaged in both epochs, included in the 
green dashed contour, is clearly the region observed with a slightly higher S/N ratio. In 
producing this catalog we also take advantage of our previous observations of the VMR-D 
region to assess if a source falls on the mapped CO-emitting cloud or not. This evaluation 
is important and interesting in this context because we do not have adequate off-source 
observations of control flelds, and the archival ^^CO observations allow us to evaluate, at 
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least to a first approximation, the possible contamination of the ON-cloud YSO population 
by OFF-cloud background/foreground objects. Operatively, we shall consider a source as 
projected on the molecular cloud, i.e. ON-cloud, if it is contained within the ^^C O contour 



l ine c orresponding to the 5 K km s ^ level in the integrated intensity map of lElia et al. 



(120071 ) which is shown superimposed on the IRAC 8 /im mosaic in Figure [2l Clearly, this is 
only a necessary condition for a source to be considered for additional constraints aimed to 
discriminate bona fide ON-cloud YSOs from galactic/extragalactic contaminants. 



3. Source extraction and photometry 

To obtain the photomet ric informatio n on the point sources detected in our mosaics. 



we adopted the DAOPHOT (jStetsoru 119871 ) photometric package. This offers a consolidated 
suite of algorithms aimed at detecting all the possible point-like objects whose signal is larger 
than a predefined threshold. The photometry was derived on our final (two epochs combined) 
mosaics. Single epoch photometry was however also measured for some sources, in order to 
assess the reliability of our final catalog. In this section we shall describe the methods used 
for source detection and photometry, as well as the selection criteria adopted for compiling 
the VMR-D point-source catalog. Completeness and reliability is also discussed and finally 
the MIPS mosaics are reconsidered to obtain a PSF photometry to be used as a complement 
to the IRAC catalog. 



3.1. Source detection 

Due to the large spatial extent of our mosaics and to the presence of extended bright 
nebulosity, distinct, local background intensity gradients are invariably present, so that it 
is not advisable to simply scan the original image with a given sensitivity threshold. The 
output of the detection algorithms in this case would be affected by a spatially variable 
sensitivity, depending on the different levels of the local background, a problem that is made 
even worse when strong gradients are present. 

To ameliorate this difficulty we first evaluated the point spread function of our mosaics 
by using relatively isolated, bright, unsaturated stars, and then performed PSF point-source 
subtractions to remove, to first order, the point-like objects. In this way we obtained images 
practically devoid of strong point sources, even if some residuals are clearly present. We 
used these images to estimate the background level across the whole field, a step that was 
done in each observed band by evaluating the sky value around each pixel in a square box 
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of 51-pixel in size. This was done by using the DAOPHOT-SKY procedure that estimates 
the modal value after clipping the outliers. Note that the relatively large size adopted for 
the scanning box is justified by the need to include a sufficient number of sky pixels even 
in crowded images, as is the case particularly for the first two IRAC bands centered at A = 
3.6 /im and A=4.5 /xm, respectively. 

Once the background in each band was obtained and subsequently subtracted from the 
original mosaics, we were left with four background-subtracted images in which the large 
scale gradients were strongly smoothed while the small scale noise was preserved. At this 
point, the compatibility of the noise with poissonian statistics is recovered by adding the 
original mean sky value back onto the background subtracted images that were then used in 
the photometry extractions. 

Source detection was then performed by scanning these final images with a star finding 
procedure which is based on a convolution with a wavelet-like filter whose spatial extent has 
been tuned to match the size of the point-like objects in the image. The signal produced 
by this convolution is then analyzed at different detection thresholds that, in our case, were 
taken at 5 a and 10 a, respectively. 

3.2. Photometry 

Once we had the list of detections from the star finding algorithm, we used aperture 
photometry to obtain a first evaluation of the source magnitudes. Although this is a relatively 
inaccurate method when crowding or background gradients are present in the image, it is 
nevertheless useful for a first estimate of source brightness, an information that was used 
later as a starting point for more accurate photometry. 

The final magnitude values were indeed obtained by means of an iterative PSF-fitting 
procedure that, in the implementation adopted by DAOPHOT, is based on an estimate of 
the instrumental PSF obtained from stellar images appearing relatively isolated in the frame. 
However, because the PSF is undersampled by the IRAC cameras, the corresponding point 
response function (PRF) that we obtain from our mosaics cannot be completely accurately 
determined; this produces a corresponding uncertainty in the subsequent fit to the stellar 
images. In spite of this difficulty, the photometry resulting from this method can be con- 
sidered more accurate than the aperture photometry because the differences in fitting the 
PRF to the stellar images remain approximately balanced between positive and negative 
residuals. In this way the total fiux is substantially preserved, even in presence of residuals, 
although at the cost of an increased uncertainty in the fitting parameters (see also § 13.41 
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below on reliability). The PRF we used is actually a combination of an analytic function (in 
our case a Lorentzian showed better results in minimizing the residuals) and a look-up table 
whose aim is to compensate for systematic differences between the pure analytic function 
and the actual stellar images. An example of the residuals remaining after subtraction of the 
modeled sources from the IRAC mosaics is shown in Figure [3l which illustrates this point 
for the two IRAC mosaics obtained at 3.6 /xm (upper panels) and the 8.0 /im (lower panels), 
respectively. 

Note, however, that the knowledge of the PRF also allowed us to take into account, 
and then correct for, the contribution to the flux of a source from wings in the brightness 
distribution of its neighboring sources. This advantage is fully exploited in crowded regions 
because the photometric values are obtained only after convergence of an iterative procedure 
that fits simultaneously groups of nearby stars instead of single objects. 

3.3. Point source selection 

We selected the sources to be included in our catalog by applying to the photometric lists 
obtained from the long exposure mosaics a sequence of logical filters based on a reasonable 
use of the quality parameters associated with the photometric procedures. 

3.3.1. Automatically selected sources 

Our first choice was to consider as "good photometry sources" only those detected above 
the 5 a threshold of the background noise. To easily distinguish faint and bright sources in 
the catalog, we also flagged those detected above the 10 a threshold. 

The DAOPHOT procedure for obtaining the point source photometry, besides evaluat- 
ing the brightness, also computes two additional quantities: CHI and SHA, related to the 
goodness of the fit and to the source sharpness, respectively. Here CHI is evaluated as the 
ratio between the observed scatter in the fitting residuals and that expected on the basis of 
the pure image noise, and is a measure of the departure of the source image from the model 
PRF. The sharpness parameter SHA, instead, is computed as a difference between the square 
of the size, in pixels, of a given source and the same quantity computed for the fitting PRF 
so that it can be used, within certain limits, as an indicator of the source angular extent. 

At first we constrained these two parameters to be CHI < 3 and liS'/Mj < 0.7 to 
represent bona fide point-like sources; we subsequently applied a further criterion to the 
flux uncertainty to select only sources with a flux relative error 6F/F < 0.5. As a third 
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step we also considered the number of iterations needed for the PRF fitting algorithm to 
converge on a given source: to avoid the worst cases, mainly corresponding to faint sources 
in crowded fields, we limited selections to those sources requiring fewer than eight iterations. 
Finally, because the PRF fitting procedure also implies a recentering of the source position, 
we filtered out all those sources in which the corresponding best fit centroid was displaced 
from the original position as evaluated in the previous detection phase by more than 1". 

3.3.2. Manually dropped sources 

Despite all our efforts to select genuine point sources, the photometric lists still included 
many detections that, on visual inspection, are actually artifacts in the vicinity of the bright- 
est sources 0. This problem affects in particular the mosaics at A = 3.6 /im and A = 4.5 /xm, 
and is related to the so called muxbleed defect. These artifacts are recognized because they 
appear as a sequence of equally spaced pseudo-sources aligned along the same row containing 
the true bright source. Because their appearance is very much like that of a train of point 
sources, these pseudo-sources escaped all our other selection filters so that it was necessary 
to eliminate them by eye. In this way more than 1000 false sources were discarded. 

3.3.3. Additional sources and short exposures 

Besides these false sources the previous selection criteria, applied to the long exposure 
photometry, discarded an additional, but relatively small, number of sources (~ 1000). These 
are mainly located in regions affected by crowding or saturation, the latter case occurring 
especially in the neighborhood of the brightest objects. To compensate for this we also 
considered the short exposure mosaics, obtaining the corresponding photometry with the 
same procedure adopted for the long exposures. The short exposure photometry was used 
when, in assembling the catalog, the long exposure photometry of a given source was rejected 
due to one of the adopted selection criteria, as opposed to pixel replacement in the mosaic 
or a hard cutoff between the short and long frames. In these cases we further checked to see 
if alternative, short exposure photometry existed that met all the catalog constraints. If it 
did, we included the source in the catalog; in this way we recovered a total of 276 sources. 

We also reexamined an additional 77 rejected sources, most of which (68 cases) are 
associated with 24 /xm MIPS sources, that escaped IRAC detection because they are either 



^See the IRAC Data Handbook, ver. 3.0, available at 



http://ssc.spitzer.caltech.edu/irac/dh/ 
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so bright they suffered saturation effects or were embedded in a bright nebulosity. For most 
of these sources we used the short exposure mosaics to obtain the photometry even though it 
is more uncertain. These objects were included in the catalog to produce a list as complete 
as possible, even though they have been flagged as "bad photometry" and are not considered 
further in the analysis of the YSO population. However the reader can use these flux densities 
with their uncetainties, but at the same time should be aware that the corresponding quality 
flags do not satisfy all the constraints previously discussed for the automatic source selection. 

In conclusion, we obtained four lists of sources, one for each IRAC band, that have been 
spatially delimited to include only objects lying in the sky region surveyed in all the four 
IRAC bands as delimited by the yellow contour line in Figure [1] 

3.4. Sensitivity, Completeness and Reliability 

An upper limit to the completeness of our catalog is clearly set by the sensitivity detec- 
tion threshold adopted in the detection phase that, we recall, was set to 5 a. This threshold 
corresponds to different limiting flux densities because the background level changes with 
the wavelength. We evaluated the limiting flux densities in each spectral band by calculating 
the average value of the pixel fluctuation in the background images obtained in the course of 
the photometric procedure. The same limiting flux densities are independently derived from 
the magnitude plots of Figure HJ showing the distribution of the catalog point sources as a 
function of their magnitude in the four IRAC bands. From these plots it is apparent that 
the logarithm of the number of sources increases proportionally to the magnitude, until the 
linear behaviour is lost and a fall in the number of sources is clearly seen. The vertical line 
drawn in each panel denotes the largest magnitude at which the catalog can be considered 
as complete, a value that is clearly lower than the limiting magnitude. These values are 
summarized in Table [1] along with their corresponding flux densities. 

As far as the reliability is concerned, we can exploit the fact that the VMR-D fleld was 
observed twice, allowing a check on the source detections. The criterion we adopt for this 
evaluation is to consider as reliable detections in a given IRAC band two possible cases: i) 
sources detected in both observing epochs, and ii) sources detected only once but showing 
at least another detection in an adjacent band. On the other hand we consider as less 
reliable detections case iii) when sources were detected only once in a given band and with 
no detection in the adjacent bands. 

To ensure that the sources we consider in estimating the reliability were observed twice, 
we restricted this analysis to the area delimited by the overlap of the IRAC mosaics from 
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both epochs (dashed green hne in Fig. [T]). Applying to these sources the previous criteria, 
we obtain the distributions shown in Figure H] as a dotted (rehable) and dashed (unrehable) 
hne. From these plots we estimate the reliability as a function of the magnitude by ratioing 
the number of sources satisfying the criteria i) and ii) with the totality of the catalog sources 
i)+ii)+iii) in a given band and magnitude interval. In this way we derive a reliability estimate 
that is always larger than 99% in the whole magnitude range. 

All the sources observed in both epochs show flux densities that are compatible within 
their uncertainities; the very few sources (aroun d 0.2% ) showing a significant photometric 
variability have been discussed in iGiannini et al.l ( 120091 ). 



3.5. MIPS photometry 

As mentioned in the introduction, we also reconsidered the MIPS observations at 24 /xm 
and 70 fim of the VMR-D region covering a slightly larger field of approxim ately 1.8 square 



degree s (see Fig. [T]). These observations have already been presented by iGiannini et al. 



( 120071 ) who reported the aperture photometry for 849 and 61 sources in the 24 /xm and 



70 fim band, respectively. 

However, because aperture photometry tends to be less accurate in the p resence of large 



background gradients such as are present in the MIPS VMR-D mosaics, in IGiannini et al. 



( 120071 ) the MIPS sources considered were limited to those in regions where reliable aperture 
photometry was possible. Because PSF fitting photometry is a more accurate tool, we re- 
analyzed the regions in the two MIPS mosaics of VMR-D to obtain a revised list of sources 
with the corresponding PSF photometry. Our aim was to obtain more accurate SEDs and 
then a more reliable determination of the spectral slopes characterizing the YSO population. 

The analysis of the MIPS mosaics was done with the same procedures described for the 
IRAC images. First, we minimized the effects of the background variations by obtaining a 
rough subtraction of the point sources from the original image. After this first step eliminated 
most of the point source signal, the mosaics were scanned with a square box of 21 pixels on 
a side to produce an estimated background image. This in turn was subtracted from the 
original mosaic producing a flattening of the large scale gradients, thus allowing a subsequent 
point source detection across the field with a much more uniform sensitivity. By using a 
5 a detection criterion we obtained two preliminary lists of sources at 24 fim and 70 fim, 
respectively, constrained by selection criteria similar to those adopted for the catalog of the 
IRAC sources (see § I3.3.ip . 



A comparison between the aperture and PSF photometry methods shows that at 24 fim 
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the aperture photometry tends to systematically overestimate the flux by up to ~50% with 
respect to the PSF photometry. Similarly, the flux differences at 70 fim are of the same 
order (< 50%), but they appear more erratic and without a systematic trend. 



In the end, we obtained two more complete lists of sources measured with PSF photom- 
etry, containing 1347 and 63 sources at 24 fim and 70 yum respectively, and we used them to 
supplement the IRAC catalog. Note that among the 24 /xm sources we include for complete- 
ness six saturated obje cts, corresponding to IRS 16, 17, 18, 19, 20, and 21 (nomenclature 
from lLiseau et al.lll992l ). In the MIPS catalog these correspond to the entries ID=146, 300, 
813, 814, 991, 1018, respectively, that had been flagged with negative 24 /xm flux densities, 
whose absolute value corresponds to the saturation limit (4 Jy). 

We adopted angular separation as the criterion for associating the 24 fim and 70 fim 
sources, requiring that the corresponding centroids be within 20"; in this way we produced a 
unique, merged set of sources with MIPS photometry. The merged MIPS catalog is available 
in electronic form (VMR-D_MIPS_catalog.txt), and it also contains information on the possi- 
ble associations of the millimeter sources as reported in previous studies. These associations 
are based on positional coincidences of the MIPS source centroid being within the beamsize 
of the mm observations. A short excerpt from the MIPS catalog is printed in Table [2j The 
histogram of the magnitudes is presented in Figure [5 and it shows that the completeness flux 
densities are essentially the same as those quoted in lGiannini et al.l ( 120071 ). namely ~ 2 mJy 
at 24 fim and ~ 300 mJy at 70 fim, respectively. Some statistics of the MIPS sources are 
presented in Table [3l 

Given this new M IPS photometry, a re-examination of the color-magnitude diagrams in 
Giannini et al.l ( 120071 . their Fig. 7 and 8) is in order. In Figure [6] we present the resulting K 
vs. K-[24] and [24] vs. [24]- [70] diagrams while in Table HJ 2nd and 3rd columns, we report the 
subdivision in the different classes based on the spectral indexes derived by using the 2MASS 
K-band magnitudes as a c omplement of the pre sent MIPS photometry. This classiflcation 
scheme is the same used in iGiannini et al.l ( 120071 ) and is based on the c orrespondenc e of th e 
colour K-[24] with the spectral slope of the different classes adopted by lRebuU et al.l ( l2007l ). 



A comparison of the present results with those in lGiannini et al.l (120071 ). based on aper- 
ture photometry, is possible. We renormalized the values in their Table 5 to reflect the 
changes in the total number of YSOs; because the measurement of the flux density now has 
changed by up to 50%, the numbers of sources in each spectral class changes as well. We 
note that: 



i) the number of MIPS sources detected now in the 24 /im band is increased while at 
70 fim it remains practically the same; 
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ii) the percentage of ON-cloud Class I sources is now ~ 10%, decreased from the previously 
estimated value of ~ 30%; 

iii) the percentage of the ON-cloud flat-spectrum sources is approximately the same, 
~ 26%, as compared with the previous estimate of ~ 29%; 

iv) the percentage of ON-cloud Class II objects is ~ 62%, an increase with respect to the 
previous estimate of ~ 37%, while Class III sources are now ~ 2% with respect to the 
previous value of ~ 4%. 

Furthermore, the more accurate modeling of the point sources allows us now to obtain 
more reliable correspondences between the MIPS, IRAC, and 2MASS objects, alleviating in 
this way the problem of multiple associations. 

In the light of the significantly lower percentage we find here for Class I spectrum 
sources, we conclude that our previous suggestion for an unusual excess of young objects in 
VMR-D does not hold up based on PSF photometry, and in all future discussions of MIPS 
sources we shall refer to our current catalog. 

4. The VMR-D point source catalog 

Once photometry was obtained we proceeded to the positional association of the sources 
detected in the different spectral bands. In this phase of band merging objects detected at 
different IRAC wavelengths were considered to be associated if the relative distance of their 
centroids was less than 1". For those cases in which, when comparing two bands, more than 
one source is found within 1", we considered only the nearest source to be associated. 

The final VMR-D catalog was produced by including all of the sources detected in at 
least one IRAC band which fulfill all the constraints reported in § 13.31 An exception to 
this rule are the 77 sources, mainly corresponding to very bright objects, whose addition we 
justified in § 13.3.31 To obtain a coherent list of sources, we restricted the final catalog to the 
field observed in all the four IRAC bands, namely the region delimited by the yellow contour 
line in Figure [H 

Obviously, by construction such a catalog contains many sources lacking good photome- 
try in some IRAC bands; only 8796 sources out of the full 170,299 have accurate photometry 
in all four IRAC bands. However, because for some purposes it is useful to have flux den- 
sisites in all bands, we choose to give a flux density value in any case. This was done by 
reconsidering less accurate photometry previously discarded in the source extraction pro- 
cess (see §3.31) or, in absence of this photometry, obtaining 5 a upper limits for the source 
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position. Consequently those missing flux densities in the catalog were recovered from the 
values that had been previously rejected due to the imposed constraints when the detection 
fall within a radius of 1" from the source position. Note, however, that all these band filling 
flux densities have been flagged as "bad photometry" in the flnal lists, in the same way 
the 5 (J upper limits are reported in the catalog for those cases without a speciflc detection 
within 1". In addition, to allow a more accurate determination of the spectral slopes and 
a more secure classification of the YSOs, the final catalog also considers the spatial coinci- 
dence between IRAC and MIPS 24 //m sources, and includes the MIPS photometry when 
the distance between centroids is within 6". A larger distance, up to 20", has been adopted 
for MIPS sources detected only at 70 jim.. These relatively large distances are close to the 
diffraction radius at 24 ^m and 70 //m, respectively, and have been adopted to include all 
the potentially associated sources. Of course, because the probability of the IRAC-MIPS 
source association decreases with the ccntroid distance, one must pay increasing attention 
when the angular distance is larger than 3"and 10"at 24/xm and 70/im, respectively. Note 
however that 70% of cases have been found within 3". With this criterion a total of 870 MIPS 
24 //m sources have been associated with IRAC sources; out of these, 689 are associated to 
a single IRAC source, 153 are double associations, 27 are triple, and one is associated with 
four IRAC sources. To allow the reader to evaluate these associations we also report in the 
catalog the distance of the centroids for each association. 

The pubhshed version of the catalog, available in electronic form (VMR-DJRAC_catalog.txt), 
summarizes the photometry arranging the information it in the following columns: 

Column (1): Source ID. Sources are sorted in ascending order of R.A. 
Columns(2-3): Equatorial coordinates (J2000). 
Columns(4-5): Galactic coordinates. 

Columns(6-7): IRAC fiux at 3.6 /xm and its uncertainty. When the fiux reported is an 
upper limit, the uncertainty has been put to zero. 

Column(8): Quality flag, assigned as follows: the "good photometry" data are flagged 
with "H" (High) if they refer to a more than 10 cr detection, and with "L" (Low) if the 
source is revealed between 5 and 10 a, respectively; in addition to this, a flag "2" or "1" 
indicates whether the source lies in the sky area covered in both observing epochs, or 
only in one of them, respectively. The "bad photometry" data (see above) are flagged 
with "b" , and upper limits with "u" . 

Columns (9- 17): Same as for Columns 6-8, but for 4.5 //m, 5.8 //m, and 8.0 //m band, 
respectively. 
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Column(18): MIPS catalog ID number of the associated MIPS source: this is reported 
if its centroid is within a distance of 6" from the IRAC source. MIPS sources detected 
only at 70 fim are reported when the centroid is within 20". Notice that, due to the 
different beam sizes, the association of more than one IRAC source with the same 
MIPS source is allowed. 

Column(19): Distance in arcsec of the associated MIPS source centroid. 

Columns(20-22): MIPS flux density at 24 /im, flux uncertainty, and quality flag similar 
to that used for IRAC flux densities but with good photometry simply flagged "G". 
The symbol NaN reported for the flux density denotes a position falling outside the 
MIPS mosaic. 

Column(23-25): The same as for column 20-22 but for the MIPS 70 /xm band. Note 
that for MIPS sources showing only the 70 /im detection (10 cases) the association is 
reported when the IRAC sources (51 cases) are within 20" from the MIPS source. 

Columns(26-27): A flag indicating the degree of associat i on wi th ^^CO(l-O) and ^^C0(2- 



1) line emission, respectively, as mapped in lElia et al.l (120071 ). "C" indicates that the 
IRAC source lies on a map pixel assigned to one or more gas clumps; "A" indicates 
that the IRAC source is not associated to a clump, but lies inside the integrated inten- 
sity contour corresponding to the 5 K km s~^, adopted here as a criterion for spatial 
association with the gas; "O" indicates a source lying inside the area observed in the 
given line, but not associated with the gas according to the two previous criteria; fi- 
nally, "N" means that the IRAC source is outside the sky area covered by observations 
in the given line. 



Colum(28): name of the 1.2 mm dust core, taken from the lists of iMassi et al.l (j2007| ) 



and reported here when the IRAC source falls within the quoted core radius. 

Examples of the catalog entries are given in Table |5] and some statistics of the catalog 
content is summarized in Table [6] and Table [71 



5. Discussion 

YSOs in SF regions can be identified by the presence of an infrared excess which is 
produced by extended circumstellar disks/envelopes surrounding these objects, and that 
participate in their pre-main sequence (PMS) evolution toward the main sequence. For 
classification purposes it is customary to use the log XFx vs log A plot to derive the spectral 
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slope, a, that is generally evaluated in the 2 /xm < A < 20 fim range. This parameter is 
taken to be dependent on the evolu tionary phase of the emitting source and its value is 



gener ally used to define four classes (iLada 1 119871 : lAndre &: Montmerle I Il994j : iGreene et al. 



19941 ): Class I for a > 0.3, "fiat-spectrum" for -0.3 <a< 0.3, Class II for -1.6 <a< -0.3, 
and Class III for a < —1.6. A further class, the Class 0, has been also introduced to include 
sources t hat, while escaping any detection at A < 10 /im, show relatively strong submillimeter 
emission (lAndre et al.lll993l ) and should then correspond to very early and cold objects. This 
relatively simple classification scheme has been h owever criticized because of its sensitivit y 
to inclination effects to the line of sight (see, e.g., iRobitaille et al.l 120071 : ICrapsi et al.ll2008l ). 

Alternative classifi cation schemes based on the bo lometric temperature ( jMyers &: Ladd 



1993; Chen et al. 1995) or the bolometric luminosity ( Young fc Evanslbood ) have also been 
proposed in an attempt to minimize the influence of the system inclination to the line of 
sight. However these alternative schemes have themselves been found to be sensitive to the 
system geometry, so that we have decided to adopt the most popular classiflcation scheme 
based on the spectral slope. This choice also allows us a more direct comparison with a series 
of studies carried out with Spitzer on Galactic star forming regions i n the framework o f the 
Legacy project From Molecular Cores to Planet-forming Disks (c2d. lEvans et al.l 120031 ). 



The reader interested in a m ore detailed discussion on the effects of adopting other clas- 
siflcation schemes should refer to lEvans II et al.l (l2009l . hereinafter E09), who also discussed 
the role of the extinction corrections. 



5.1. Color analysis, and YSO selection 

In the following analysis we focus our attention on the 8796, out of a total of 170,299 
catalog sources, for which good photometry was obtained in all the four IRAC bands; i.e. 
the working sample already deflned in ^ In addition to this database we also consider, 
when available, the /T-band 2MASS photometry and the 24 fim MIPS photometry to allow 
a more accurate estimate of the spectral index in the 2-24 fim spectral interval. 

Our flrst step was the determination of this spectral index a from the linear flt of the 
observed log XFx vs log A points for all the working sample sources. The distributions of the 
a values obtained are shown in the two panels of Figure [3, representing the spectral indexes 
for ON-cloud (upper) and OFF-cloud (lower) sources, respectively. In the two panels the 
histograms are drawn separately for the objects in the working sample (solid line) and for the 
fraction of the candidate YSOs (dashed line) that have been subsequently selected according 
to the criteria we shall discuss in the next section. 
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We note that the number of objects falhng within the ON- and OFF-cloud region are 
quite similar (4484 against 4312). Such a fact derives from the similarity of both areas 
(0.57 against 0.60 squared deg) and indicates that the overall population of both regions 
is dominated by photospheres and background/foreground objects, while the YSO popula- 



Evans II et al. 


2009; 


Rebull et al. 


2010 



Guieu et al.l 120091 ). As a consequence, stringent criteria have to be exploited to select the 



candidate YSOs. 



5.1.1. Source selection 



T he ab ility of the IRAC colors to ide ntify YSOs has been st udied by iHarvey et al. 



(120061, 120071 ) in the Serpens cloud and by iGutermuth et al.l (120091 ) in a large number of 
YSO clusters, and we us e parts of both a pproa ches here, augmented by criteria identifying 
AGB contaminants from iMarengo et al.l (120081 ). After considering the locus of the normal 
reddened stars and the non-negligible extragalactic contamination, these authors suggested 
two statistical criteria to select candidate YSOs in the [8] vs ([4.5]-[8]) col-mag diagram, 
namely: i) ([4.5]-[8]) > 0.5 and ii) [8] < 14- ([4. 5] -[8]). This first criterion excludes the 
stellar component while the second, designed to statistically cut the 95% of the extragalactic 
objects, has been tuned by studying the extragalactic component observed with IRAC in 
the Elais Nl SWIRE region. 

Applying these criteria to the VMR-D cl oud we select 595 YSO ca ndidates ON cloud 
and 244 OFF-cloud, respectively. Following iHarvey et al.l ( l2006l . 120071 ). to include other 
potential colder candidates that could have been escaped in this diagram, we also considered 
further diagrams based on longer wavelength photometry, i.e. the [24] vs ([8.0]-[24]) and 
[24] vs ([4.5]-[8.0]) diagrams. Adopting the appropriate selection criteria we included 25 
more sources ON and 17 OFF-cloud. The YSO candidates selected are shown in Figure H] 
(ON-cloud) and (OFF-cloud) as open circles and grey crosses, representing objects showing 
Class I/flat and Class II/III SEDs, respectively. 



A lso in dicated on these plots are the original selection criteria from the iHarvey et al. 



(l2006l . 120071 ) study of the Serpens cloud. However, note that the Ser pens cloud (260 pc) 
is much closer than the VMR-D cloud (700 pc, see iLiseau et al.lll992l ). To reflect this, we 
have also indicated in the figure the region into which a candidate YSO in Serpens would 
move if displaced to the VMR-D distance. The region between the two inclined parallel 
lines can then be regarded as potentially containing both candidate YSOs and extragalactic 
contaminants meaning that objects in this region should be more accurately studied to 
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ascertain their nature. Here we conservatively consider only sources delimited by the original 
criteria represented by the dashed lines in Figure [8] and [91 

In the same figure the two color ([3.6]-[4.5]) vs ([5.8]-[8.0]) diagram is also shown in the 
upper right panel with the same meaning of the symbols, and with a square box that i s 
drawn to show the approximate locus of the Class II sources according to lAUen et al.l ( l2004j ). 
The irregular area del imited by the dot-das hed line is the locus of the AGB stars in this 
diagram, according to iMarengo et al.l (120081 ) . For comparison, in a separate Figure [HI we 
show the same diagrams drawn for the OFF-cloud sources. 

With the aim to obtain a list of candidate YSOs as clean as possible f rom contaminants 
we als o adopted further criteria recently discussed and summarized by iGutermuth et al. 
( 120091 ). These have been derived for identifying as likely contaminants those sources satisfy- 
ing a series of constraints flagging star-forming galaxies and AGN, as well as shock excited 
knots that can mimic point sources in star forming environments. These are based on color- 
color diagrams as [4.5]-[5.8] vs [5.8]-[8.0], [3.6]-[5.8] vs [4.5]-[8.0], and [3.6]-[4.5] vs [4.5]-[5.8], 
as well as on the color-magnitude [4.5] vs [4.5]- [8.0] diagram. 

In Figure [To] we present these diagrams (displayed for simplicity only ON-cloud) for the 
sources emerging from the previous selecti on step that are reporte d with crosses and circles 
in Figure [H Considering the constraints in IGutermuth et al.l (120091 ) we flagged 64 sources, 42 
ON (shown as open squares) and 22 OFF-cloud, respectively, that have then been eliminated 
from our list of candidate YSOs. 



As a further check we also applied to this list the IRAC based criteria of [Marengo et al. 



( 12008[ ) for flnding AGB contaminants, identifying in this way 115 sources ON-cloud (112 
showing a Class III SED) and 65 OFF-cloud (62 showing a Class III SED). Also these 
sources have been eliminated from the flnal list. 

At the end of this process we remain with 637 IRAC selected candidate YSO sources, 
463 ON and 174 OFF-cloud that we assign to a speciflc class simply considering their spectral 
slope as derived considering all the available data in the 2-24 fim spectral interval. These 
sources have been used to draw the dash-dotted histogram in Figure [3 The corresponding 
distribution in Classes of these candidate YSOs is shown in 4th and 5th column of Table [H 



5.1.2. Comparison with other Star Forming Regions 



The analysis presented above allowed us to classify the VMR-D young population and 
it is interesting to compare our results to those obtained by Spitzer for other SF regions. 
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These have been summarized and analyzed in a consistent way by lE09l and we adopted it 
for the VMR-D comparison. Ho wever, because our method for selecting the YSOs is not 
exactly the same as that of IEOQI . before to compare our results it is desirable to evaluate 
the e ffects produced by applying our pipeline at least to one of the five clouds considered by 
EOQI . The Serpens cloud is probably the most appropriate for this check: it is projected near 
the Galactic plane and has bee n scrutinized with a ccurate spectroscopic follow-up aimed 
to identify AGB contaminants (jOliveira et al.ll2009l ). For this cloud we do not attempted 
to select candidate YSOs out of the full catalog, dehvered by the c2d legacy project and 
made available at the Spitzer Science Center, because our selection criteria are similar to 
those used to obtain the c2d YSO list. The difference in our approach is mainly in the 
identification of the contaminants so that we used the Serpens YSO catalo g to analyze with 



our pr ocedure the 262 sources listed. Applying our pipeline, based on the iGutermuth et al. 



(I2OO9I ) criteria, to this list we exclude a total of 14 s ource s. None of these is in the list of 
the 78 candidate YSOs scrutinized by lOliveira et al.l ( 120091 ) and their distribution in classes 
is: two sources in Class I, two are flat-spectrum, eight in Class II, and two in Class III. 
Despite we exclude these 14 sources out of the 262 in th e Serpens YSO catalog, we note that 
after their subtraction from the numbers quoted in lE09l and reported in Table El the relative 
distribution in classes remains practically the same and also the star formation rate is not 
significantly affected. In the limit of this in dicat ion we can compare in Table [8] the results 
obtained for VMR-D with those provided by lE09l . n oting that in the nurn bers quoted for the 
YSOs in Serpens the sources identified as AGB by lOliveira et al.l (120091 ) have been already 
subtracted. 



In this table the rows 1-3 consider the sources selected by means of the iHarvey et al. 



(I2OO7I ) criteria and shown in Fig. |8] and IH The subsequent rows 4-6 report the results 
obtained after filtering the sources in the first thre e rows with more stringent criteria for 



elimi nating extragalactic and galactic contaminants ( iGutermuth et al.ll2009l : iMarengo et al. 



20081 ) as illustrated in Fig. [101 Further three rows, namely 7-9, are used to present the results 
obtained after apply ing to the first three rows the correction for contaminants suggested by 
Oliveira et al.l (|2009[ ) on the basis of their spectr oscopic investigations of the Serpens YSO 
population selected with the iHarvey et al.l ( 120071 ) criteria. 



Indeed, because Serpens is a SF region quite close to Galactic plane (/ ~ 5°), it is 
reasonable to assume that its level of contamination of AGB and giant stars is comparable 
to that of VMR-D. The sample of YSOs in Serpens is contaminated by background giants 
to a level of about 25% so that, by applying the same correction to our VMR-D ON and 
OFF-cloud sample, we obtain the statistics given in Table [HI lines 7-9. Note that in this 
case we do not report the distribution in classes because the 25% contamination is referred 
to the total number of YSOs. 
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In the subsequent lines 1 0-14 we also report for comparison the values of the corre- 



sponding quantities quoted by IE09I for Chall, Lupus, Perseus, Serpens, and Ophiucus. In 
each line we give the investigated area, the adopted cloud distance, the total number of YSOs 
found, the surface density per solid angle and per square parsec, the number of objects of 
type Class I, flat-spectrum. Class II, and Class III, respectively, and finally the derived star 
formation rate. 

First of all we recall that, differently from these other SF regions, VMR-D is located 
in the Galactic plane and it is also from 2 to 5 times farther away. As already mentioned, 
the latter circumstance imposes a larger value for the limiting luminosity (typically 0.1 Lq 
instead of, e.g., 0.004 Lq for the Taurus) resulting in significant contamination from stars; 
conversely, it offers the opportunity to study the SF process where the mass reservoir is the 
largest available. 

As a rough first approximation, we consider the OFF-cloud statistics as representative 
of the average population of the Galactic plane and, as such, it could be subtracted from 
the ON-cloud statistics after weighting for the different areas, in this way estimating the 
uncontaminated young population of VMR-D. By doing this we obtain the lines 3, 6, and 9, 
of Table [HI with the caution that the OFF-cloud region could be active in SF as suggested 
by the the histograms in Figure U\ albeit at a lower level. Because of this and of the larger 
luminosity limit, the numbers reported here for the YSOs should actually be considered 
lower limits. However, taken at their face value, they suggest an excess of Class III (~22%) 
and a deficit of Class II objects (~48%) with respect to the other SF regions. The cause 
of the apparent excess of Class III might be either evolution - the Ve la cl oud is older than 



these others - or due to uncorrected contamination of AGB stars (see lE09l ). 



Considering the statistics corrected for extragalactic and galactic contaminants, reported 
in the 4th and 6th lines, make the distribution in classes of the YSOs in VMR-D resemble 
more that seen in Serpens and Ophiuchus rather than in other SF regions. Differences 
actually exist in the Class II/III relative ratios, but they probably stem from an incomplete 
correction since VMR-D is more embedded in the plane than is Serpens. If we try to account 
for such discrepancies, we can reasonably rule out evolutionary effects since the relatively 
large fraction of Class III sources in VMR-D does not reflect a particularly small fraction 
of Class I and fiat-spectrum sources that instead appear relatively numerous. Moreover, if 
we calculate the SF rate per square parsec, by assu ming a mean mass of 1 Mq, a fraction 



of binary of 0.5 and a period of 2 Myr for SF (jE09l ). we obtain for the VMR-D SF rate ~ 
2.5-4.1 M0Myr~^pc~^, a range of values that is more consistent with that of the younger 
clouds like Serpens, Ophiucus, and Perseus, than it is for the smaller corresponding values 
obtained for Lupus and Chamaleon. 
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It is also noteworthy that adopting the corrections suggested by lOhveira et aLl (j2009[ ) 
we obtain (Tab. [8|, hnes 6-8) very similar results to those obtained after considering the set 
of criteria discussed in § 15.1.11 to exclude likely contaminants (lines 4-6). This adds some 
confidence in the present results. 

Given these observational indications it seems plausible to consider the VMR-D age as 
being in the ra nge 1-2 M y r, in good agreement with the dynamical time estimate of 1.5 
Myr derived by 



flPalla fc Stabler 



_lia et al.l (120071 ) fo r VMR-D and comp arable to the estimates for Perseus 
2000h and Serpens (IDjupvik et al.ll2006l ). Within the previous assumptions, 



and given the solid angle subtended by our ON-cloud region, it is easy to estimate that ap- 
proximately 200-300 Mq Myr~^ are converted into s tars in the inves tigated region. Because 
the cloud mass has been estimated at 1.5 x 10^ Mq (lElia et al.ll2007l ). we can also derive the 
depletion time for the cloud, that is tdepi ~ 40-70 Myr, as well as the SF efficiency given 
by MYso/(MYso+Mc;oMd)~ 0.02-0.03. However, the uncertainty associated to the estimate 
of the cloud mass is at least AM/M ~ 50% so that these quantities are correspondingly 
uncertain. 



5.2. IRAC vs. MlPS-based classification and final lists of YSO 

We have seen that the IRAC data are an essential, additional ingredient of YSOs clas- 
sifi cation, and by i ncludi ng them we have been able to reconcile (as indeed was anticipated 



by iGiannini et al.l (120071 )) the estimated excess of Class I objects previously derived using 



only MIPS aperture photometry. 

In this section we want to briefly discuss the young population census, as obtained by 
means of the present IRAC photometry complemented with 2MASS and MIPS data, with 
the results that can be derived using only 2MASS and MIPS photometry. The aim is to 
assess the ability of the MIPS photometry to correctly classify YSOs when no IRAC data 
are available, as is the case for the MIPS regions without overlap with the IRAC mosaics 
(see Fig. [ID . 

The classification method used in these two cases is different: in the first case the 
2MASS/IRAC/MIPS photometric data are used to evaluate the spectral slope a and to 
separate candidate YSOs from contaminants, while in the sec ond case only 2MASS /MIPS 



data are used to exploit the vs. Ks-[24] diagram as in iGiannini et al.l (j2007| ). The 
first method is clearly the most appropriate involving more photometric information (up 
to six wavelengths); however, being the second method the alternative choice when MIPS 
observations are the only available, it is interesting to check if differences arise in selecting 
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the YSO population with these two methods. The two hsts of candidate YSOs obtained 
with these two approaches are reported in the corresponding electronic catalogs: VMR- 
D_IRAC_YSO.txt and VMR-D_MIPS_YSO.txt. Typical entries for both are reported in 
Table M and [101 respectively. 

In discussing the differences emerging between these two approaches we shall limit our- 
selves to consider only the corresponding lists of Class I objects; similar comparisons for the 
other Classes (flat, II or III) follow the same line of thinking. Before proceeding, we reiterate 
that the census of the different Classes in this paper is based on the working sample, i.e. is 
limited to those sources with valid fluxes in all four IRAC bands. This clearly excludes the 
sources with an upper limit, a large flux uncertainty, or even a not strictly point-like shape 
in at least one IRAC band. These objects, although not considered here, are included in the 
general catalog and could represent interesting targets for future investigations. 

The 49 sources selected as Class I objects by scrutinizing the IRAC catalog working 
sample are listed in Table [TTl Along with the source ID, coordinates and spectral index, we 
also report the presence of a counterpart in other bands: K<j-2MASS in column 5 and MIPS 
bands in columns 6 to 8. This information is needed to compare the two methods since the 
MIPS classification is based on K<,-2MASS and 24 fim magnitudes, and only IRAC sources 
associated with both these ancillary data are subject to both IRAC and MIPS classification. 
By looking at Table [TT] (column 9), we see that only five sources out of the 49 entries were 
classified with MIPS and all of these are Class I (flag I). While this is the expected result 
when the four IRAC points are reasonably well aligned with the two extreme spectral points 
at Ks and 24 /im, it is however noteworthy that in the relatively few cases in which both 
methods can be applied they give the same classification. In columns 10 and 11 we also 
indicate whether the source is associated with gas and dust emission. In this respect, the 
genuine nature of the selected Class I sources is testified by the fact that 44 out of the 49 
sources are coincident with gas clumps or peaks. The remaining 5 sources are located in 
regions not covered by the ^^CO map (flag N). Finally, the last column in dicates whether 



we fi nd that the source is also IRAC variable on a timescale of 4.5 months ( IGiannini et al. 



20091 ): we identify six cases of variability. 



Table [T^ is complementary to Table [TTl It provides the list of all the 20 sources classified 
as Class I based on Ks-2MASS and 24 yum photometry. Besides the Ks-24 index (column 4) 
information is also given about the presence of a 70 fim counterpart (column 5). The 
association with an IRAC counterpart is indicated in column 6 while column 7 reports the 
quality flags of the detections in the four IRAC bands. Since this classification is determined 
by MIPS, there are only seven sources (see column 8) contained in the IRAC working sample, 
while 11 further sources are associated with IRAC objects lacking good photometry in one 
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or more bands. The remaining two sources do not have a counterpart in the IRAC catalog 
(see the notes to Table [T2l) . 

Five objects in this table are clearly in common with Table [T^ but there are two more 
cases included in the working sample and classified as fiat-spectrum sources when we derive 
the spectral slope exploiting the IRAC data. This suggests that the classification inde- 
pendently provided by K,-2MASS/MIPS and K,-2MASS/IRAC/MIPS is not completely 
equivalent even if the two methods seem to be in agreement in five out of the seven common 
cases encountered. Again, most of these sources are located in gas clumps. Note also that 
Table [12] contains many sources whose IRAC fiags do not fulfill our compelling requirements 
(only HI and H2 quality fiags). Although they are not members of the working sample, they 
are nonetheless good candidates for very young protostars. 

Finally, for the sake of completeness, we present in Table [12] a list of MIPS selected 
sources that, for various reasons, can be considered as particularly red or cold. When 
observed in both MIPS bands, these objects have a color [24]- [70] > 6.3, a condition we 
imposed to select the coldest (T < 50K) sources (see Figure [HI right panel). Other sources 
listed in Table [13] are: i) the 70 fim sources with point-like appearance but without a 24 fim 
counterpart; ii) the objects saturated at 24 fim, but already classified as IRAS Class I objects 
(Liseau et al. 1992). Note that increasing the angular distance of the association, reported 
in column 7, the uncertainty also increases so that further analysis is required to confirm a 
genuine source association. 



5.3. Clustering 



The spatial distribution of the young stars provides useful information about SF mech- 
anisms in the cloud, since YSOs do not move very far from their birthplaces during the PMS 
evolution. Much of the accumulated observational evidence suggests that stars prefer to form 
in groups, and in particular, that the two-point correlation function of youn g stars show s 
a different behaviour at small and large spatial scales, respectively (see, e.g., ISimonll 19971 ). 



More recently iGutermuth et al.l (120091 ) analyzed the spatial distribution of the YSOs in 36 
clusters by means of both the nearest neighbor and the minimum spanning tree method. A 
critical length, in the range ~ 0.1-0.5 pc, could be identified in the cumulative distribution 
of the source spacings indicating differences in the slope at small and large spatial scales. 
We use here the function w{6) to quantify the degree of the YSOs clustering in the VMR-D 
cloud: 



w{e) 



No{e) 



1 
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where Np{6) is the average number of pairs with angular separation between 6 and 6 + A6 
around each object, and Nq{6) is the same number expected for a random distribution of 
the same sources. For SF regions, this function for PMS objects typically shows a power 

law behaviour w(9^ OC ^ but with cle^i"!'^'^ HifFprpnt p-ynnnpnts for small ex 9"! nnH larcrp 

(7 ~ 0.5) spatial scales, respectively 
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This fact has been interpreted as an indication of a critical scale that separates the regime of 
binary and multiple s t ars, a t smaller scales, from the hierarchical clustering dominating at 



larger scales. iLarson I (119951 ) also suggested that this scale could represent the Jeans length 



of the molecular cloud, although this interpretation encounters difficulties in explaining th e 



different length scales observed in Taurus, Ophiuchus and Trapezium clusters (ISimorull997l ). 



We applied the w{6) statistics to our candidate YSOs selected among the VMR-D 
working sample sources (see § 15.11) : that is, those sources satisfying the criteria for minimizing 
the presence of galactic and extragalactic contaminants illustrated in § 15.1.11 and Figures [H 
m and[TOl The spatial distribution of these sources is presented in Figure [2] where the objects 
of Class I/fiat (red circles) and Class II/III (green crosses) are shown superposed to the 8.0 
yum IRAC mosaic. 

We find substantial differences between clustering in early and late Classes as is shown 
in Figure [11] where we present the behaviour of the w{9) function for Class 1 and flat- 
spectrum (left panel), and Class 11 and III sources (right panel), respectively. It is apparent 
that in both cases the correlation function behaves as a single power law, at least in the 
investigated range of scales, but with noticeably different slopes denoting a higher degree of 
clustering of the younger classes with respect to the later evolutionary stages. This can be 
also qualitatively seen by eye from the spatial distributions in Figure [2l 

Note, however, that in the range of the investigated spatial scales, there is no break in the 
power law slope withi n each region- this is significant even though a direct comparison with 
the cases discussed in ISimonI (119971 ) for optically selected PMS stars is not straightforward 
because here we deal with IRAC selected YS Os instead of optica. lly selected PMS stars. A 
more appropriate comparison can be done with lHarvey et al.l (120071 ) who show the behavior of 
the same w{6) function computed for the distribution of the YSOs in Serpens (their Figure 
10). The slopes they found for both Class I/fiat and Class II/III sources are remarkably 
similar to those shown in Figure [TT] for VMR-D sources. A possible difference may be due to 
a trend in Serpens toward smaller scales due to a flattening of the function at scales smaller 
than 5 X 10^ AU. 

A plausible justification for the different slopes of early and late YSOs is simply that 
it is due to the different ages of the two populations. It is generally accepted that stars are 
essentially formed in clusters, so that it is naturally expected that the older objects have 
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moved from their birthplaces, producing object dispersal we measure in the shallower slope 
of the correlation function for Class II/III sources. 

Finally in Figures [T2] and [12] we show more details of the four regions around IRS16, 17, 
19, and 20, including the spatial distributions of the Class I/flat (red circles) and Class II/III 
(green crosses) located within 4'. To produce these figures, we consider an enlarged sample 
containing all the catalog sources located around the IRS objects and with good photometry 
in at least two adjacent IRAC bands. This allows us to compensate for the relative scarcity 
in these small areas of candidate YSOs selected from our working sample and then to better 
highlight the trend in the spatial distribution of the different classes. The sources reported 
here have been selected from this enlarged sample by applying the same criteria discussed 
in § 15.1.11 and illustrated in Figure El El and [101 In doing this we discarded the sources 
with upper limits and considered only those cases for which at least a bad photometry was 
obtained in the two complementary IRAC bands. 

In all cases we see that around these four IRS objects the Class I/flat sources (red circles) 
are not uniformly distributed. They appear more abundant around IRS16 and IRS17; Class 
II/III objects (green crosses) are instead more numerous and more simmetrically distributed 
in IRS16, 19, and 20. An exception is IRS17 where the Class II/III objects follow the same 
elongated distribution of the Class I/flat sources. We also note that this configuration in 
IRS 17 is suggestiv e of the presence of two separated clusters as has also been suggested by 
Massi et al.l (12006! ). Speculating on this case, we suggest that a dust disk, whose plane is 
perpendicular to the observed YSOs distribution, may be responsible for this appearance; 
this is suggested also by the central dark lane visible in all the IRAC bands. It is however 
clear that a conclusion cannot be given here and further investigations are needed. 



6. Conclusions 

This paper presents the analysis of the Spitzer-IKAC mosaics obtained in the region 
of the Vela Molecular Ridge, Cloud D. In a ~ 1.2 square degrees area, corresponding to 
the overlap of the IRAC mosaics in the four bands, we obtained the photometry that has 
been used to compile a four-band point source catalog. This has been supplemented by 
relevant information at longer wavelengths from revised photometry of the VMR-D MIPS 
mosaics and from previously published millimeter line and continuum observations of the 
same region. This catalog has been used to outline the general characteristics of the YSOs 
in this SF region. The main results are: 

1. The VMR-D IRAC point source catalog has been compiled with 170,299 entries, 8796 
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of which, constituting our working sample, have good quahty fluxes in all four IRAC 
bands. Complementary information has been also added using a revised, PSF-based, 
MIPS 24 fim and 70 /xm photometry and previous millimeter observations obtained in 
both CO lines and dust continuum. 

2. A revised catalog containing the point sources extracted from the VMR-D MIPS 24 fim 
and 70 fim mosaics has been also obtained following the same photometric procedure 
adopted for the IRAC bands. 

3. The working sample has been analyzed according to the spectral slope that was ob- 
tained by also considering both the Ks-2MASS and the MIPS 24 fj,m fluxes when 
available. 

4. We selected the candidate YSOs by means of their position in color-color and color- 
magnitude diagrams involving IRAC and MIPS photometry; the likely contaminants 
due to extragalactic and galactic objects have been also estimated and excluded from 
the final list of YSOs constituted by 637 objects. 

5. The candidate YSOs have been used to evaluate the surface density (3.3-5.5 pc~^), 
the star formation rate (2.5-4.1 M0Myr~^pc~^) and efficiency (~ 0.02-0.03), and the 
depletion time (40-70 Myr) of the VMR-D cloud. 

6. The global clustering properties of the YSOs, studied by means of the t wo point cor- 



relatio n function, have been found to be very similar to those derived by iHarvey et al. 



(120071 ) for the Serpens SF region; it clearly shows that Class I/fiat sources are more 



clustered than Class II/III sources. 

The IRAC catalog sources around IRS 16 and IRS 17 show a larger fraction of Class 
I/flat SEDs with respect to IRS19 and IRS20 where the Class II/III objects are rela- 
tively more numerous. In the IRS17 region there is a tendency for the Class II/III to 
trace the same elongated distribution of the Class I/fiat sources. 
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Fig. 1. — Three-color image of the observed field, obtained by superposition of the mosaics 
acquired in different bands. The color coding is blue, green, and red for A = 3.6, 5.8, and 
8.0 fim, respectively. A yellow contour line delimits the region imaged in all the four IRAC 
bands. The geometry of the scans, obtained in two different epochs, is such that only the 
central part of the field (dashed green line) is observed twice, while more external regions 
are covered only once. The more external red dashed lines delimit the fields observed with 
MIPS at 24 /im (upper) and 70 /im (lower), respectively. 
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Fig. 2. — The ^^CO contours from lElia et al.l (120071 ) are reported (yellow line) superposed 
to the IRAC 8.0 /im mosaic, starting at 5 K km s~^ in steps of 15 K (km s~^). The yellow 
dashed line delimits the region observed in ^^CO. The spatial distribution of the Class I/flat 
(red circles) and the Class II/III (green crosses) is also shown. The reported sour c es are 
those recognized as bona fid e YSOs by using the se lection criteria of iHarvey et al.l (|2007| ) 
and the exclusion criteria of iGutermuth et al.l (|2009[ ) discussed in § 15.1.11 and illustrated in 
Figure El El and HDl 
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Fig. 3.— Portion of VMR-D cloud delimited by a=(08'^ 47"" 45") and 5=(-43°57'54") at 
the lower left corner, and by a=(08^ 47'" 36*) and 5=(-43°56'41") in the upper right. The 
left panels show the mosaics obtained at 3.6 fim (upper) and 8.0 fim (lower), respectively. 
The right panels illustrate the residuals remaining after subtraction of the point sources that 
have been modeled with the PRF obtained by using isolated unsaturated stars in the same 
mosaics. 



Table 1. Overview of the limiting magnitudes and completeness 







IRAC 




MIPS 




3.6 


4.5 5.8 


8.0 


24 


70 


"p a 
^ lim 


1.1 X 10-^ 


1.0 X 10-^ 3.5 X 10-^ 


4.9 X 10~^ 


6.3 X 10"^ 


1.0 X 10-^ 


miim 


18.5 


18.1 16.3 


15.3 


10.1 


2.2 


F a 
comp 


4.9 X 10-^ 


4.5 X 10-^ 2.9 X 10"^ 


4.0 X 10-^ 


2.0 X 10-3 


5.0 X 10-^ 


^comp 


16.9 


16.5 14.0 


13.0 


8.9 


1.0 



^AU fluxes are given in Jy 
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Fig. 4. — The distribution of the magnitudes in the VMR-D catalog. Each panel corresponds 
to an IRAC band, as indicated. The upper curve (continuous hne) represents all the sources 
included in the catalog. The dotted line refers to sources observed in both observing epochs 
(i.e. falling within the green dashed line in Figured]). Here the sources included are either 
detected twice or even once but with a detection in an adjacent band. The lower dashed 
curve corresponds to sources considered unreliable because they are detected only once, 
without any detection in adjacent bands. From these distibutions the catalog completeness 
and the reliability are estimated (see text). The vertical dashed line is drawn to mark the 
completeness limit. 
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Fig. 5. — Flux histograms for the VMR-D point sources detected in the MIPS mosaics 
at 24 //m (left panel) and 70m //m (right panel). In the left panel, 6 saturated sources 
corresponding to IRS objects have been discarded (see text). 
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Fi g. 6. — Color-m a gnitu de diagrams for MIPS sources. A comparison with similar diagrams 



m 



Giannini et al.l ( 120071 ) shows that by adopting PSF instead of aperture photometry, a 
systematic decrease (< 50%) of the 24 /xm flux and random variations (to a lower extent) 
of the 70 /xm one has occurred. As a consequence, the number of Class I sources, estimated 
through 2MASS and 24 /im photometry (left panel), is now diminished with respect to the 
early prediction (see also § 15. 2p . Moreover, because of the correspondingly larger [24]- [70] 
color, an increased number of objects now populates the area corresponding to T < 50K, in 
the right panel. 
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Fig. 7. — The distribution of the spectral index a, as derived by fitting the source fiuxes be- 
tween 2 and 24 fim (see text). In both panels the continuous line refers to the 8796 "working 
sample" sources (4484 ON-cloud, and 4312 OFF-cloud) with good photometry in the four 
IRAC bands. After subtraction of both normal stars and extragalactic contamination (see 
text and Figures [HI [9] and [TOj) we remain with 637 sources, that are located either ON-cloud 
(dashed line, 463 cases) or OFF-cloud (dashed line, 174 cases). 
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Fig. 8. — Color-magnitude and color-color diagrams for the VMR-D point sources detected 
in the four IRAC bands and located ON-cloud (see § [2]). Sources selected applying the 



Harvey et al.l ( l2006l 120071 ) criteria are plotted as circles (spectral index a > —0.3) and grey 
crosses (a < —0.3). In the upper panels grey dots represent the 8796 "working sample" 
sources. The reddening effect, c orresponding to a visual extinction of Av=50 mag, is shown 
by the arrow (Ilndebetouwll2005[). In the col o r-magnitud e plots the dashed lines correspond to 
statistical criteria derived by lHarvey et al.l ( l2006l . 120071 ) for separating YSO candidates from 
normal stars (vertical lines) and extragalactic contamination (oblique lines) in the Serpens 
dark cloud. The dash-dotted line represents the displacement of the dashed line after scaling 
to the VMR-D distance (see text). The dotted lines show the completeness limit and the 
horizontal dashed line shows the 24 /im magnitude beyond which a source is considered 
extragalactic. In the two-color upper right panel the small s quare shows the approximate 
domain for the Class II sources according to lAUen et al.l ( 120041 ). The dot-dashed line encloses 
the region corresponding to the AGB colors, but note that a source can be considered as a 
candidate AGB only when two further constraints are satisfied in the [3 .6]- [4.5] vs [3.6]-[8.0] 
and [4.5]-[5.8] vs [5.8]-[8.0] diagrams, respectively ( jMarengo et al.ll2008l ). 
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Fig. 9.— 



The same as in Fig. |H] but for sources located OFF-cloud. 
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Fig. 10. — Diagrams use d to isolate the contamina nts among the sources of the "working 
sample" selected with the iHarvey et al.l ( l2006l . 120071 ) criteria. For simplicity only ON-cloud 
sources are displayed. In each diagram, solid lines delimit the loci where possible contami- 



nants are located (see iGutermuth et al.ll2009l ). The grey circles represent the sources shown 



in Fig. [8] as crosses and circles. The upper panels are used to isolate unresolved starforming 
galaxies falling under the continuous line. The lower left panel isolates broad-line AGN and 
the lower right is used to identify likely shock emission (upper left region) and possible cases 
of PAH contamination (lower right). Objects falling in at least one of the regions delim- 
ited by the continuous lines are identified as likely contaminants and are reported as open 
squares. 
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Fig. 11. — The two-point correlation function for all the candidate YSOs in VMR-D cloud. 
Left panel: sources classified as Class I and flat-spectrum. Right panel: Class II and III 
sources. The number of objects is reported in each panel with the slope 7 of the linear fit. 
The slope uncertainty, evaluated by bootstrap resampling, is ±0.10 and ±0.08 in the left 
and right panel, respectively. 



[47:59.9 . i 

-_f 



Oi.2 . / 5(^.4 iA_Ai6____ 



-JW-l-ir-^-S^i^— : mk --:ie:7 



Fig. 12. — Class I and flat-spectrum (red circles) and Class II and III (green crosses) IRAC 
sources found in a radius of 4' around IRS16 (left panel) and IRS17 (right panel), respectively. 
Sources shown here possess at least three IRAC flux densities of good quality (see text). 
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Fig. 13. — As in Figure [12] but for the regions around IRS19 (left panel) and IRS20 (right 
panel), respectively. 
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Table 2. VMR-D MIPS catalog: typical entries. 



A/TTTDQ TT^ 


KA(ZUUU j 


UeC[ZV\)\j ) 


Lrion 




r z4 




r 




tjlij 


At/ 




OU(l-0j 


OU(2-lj 


Uust 




1 


130.79321 


-44.01385 


263.47831 


-0.97396 


8.420e-03 


1.812e-03 


O.OOOe-fOO 




O.OOOe+00 


0.0 


N 


N 


N 






2 


130.81630 


-43.97870 


263.46093 


-0.93919 


1.535e-03 


1.135e-04 


O.OOOe+00 




O.OOOe+00 


0.0 


N 


N 


N 






135 


131.38707 


-43.71376 


263.50867 


-0.45199 


4.314e-03 


1.037e-04 


O.OOOe+00 




O.OOOe+00 


0.0 


I 


N 


N 






136 


131.38847 


-43.83051 


263.60060 


-0.52397 


1.934e-^00 


2.598e-01 


2.001e+01 




8.478e-01 


2.1 


I 


C 


C 


MMSl 





^Coordinates are given in degrees; fluxes (F) and their uncertainties (E) are given in Jy. 
'^Angular distance in arcsec between 24 /im and 70 /im centroids. 

'^I and N are used to flag a source as falling inside or outside the IRAC observed region. 
'^Symbols are explained in §[4] (Columns 26-28). 



Table 3. MIPS source statistics 



CO association 


24 /im only 


70 yUm only 


24 & 70 /im 


Total 


Inside contours 


604 


7 




651 


Outside contours 


131 


1 


2 


134 


Outside IRAC catalog 


560 


8 


10^ 


578 


Total 


1295 


16^ 


52 


1363 



'^Five sources are saturated at 24 /tm corresponding to IRS16, 17, 19, 20, 
21; five other sources have a duplicate 70 /im association. 

'^One source is saturated at 24 /im, corresponding to IRS18. 

*^Four sources are outside the 24 /xm mosaic; visual inspection of the re- 
maining shows that nine are diffuse while three appear more compact and are 
reported in Table [12] as remarkable objects. 
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Table 4. Classification based on MIPS and IRAC photometry 





MlPS-based 


IRAC-based 


Class 


ON-cloud 


OFF-cloud 


ON-cloud 


OFF-cloud 




N 


N 


N (%) 


N (%) 


Class I 


16 (10.1) 


4 (3.4) 


44 (9.5) 


5 (2.9) 


Flat spectrum 


41 (25.9) 


20 (16.9) 


61 (13.2) 


10 (5.7) 


Class II 


98 (62.0) 


84 (71.2) 


230 (49.7) 


93 (53.4) 


Class III 


3 (1.9) 


10 (8.5) 


128 (27.6) 


66 (37.9) 



'^These MlPS-based p e rcenta ges can be compared with those in 
Tab. 5 of ICiannini et al.l (120071 ) after renormalization of their values 
to the YSOs only. 
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Table 5. VMR-D IRAC Point Source Catalog: typical entries 



ID 


RA(2000) 


Dec(2000) 


Glon 


Glat 


F3.6 


E3.6 


Q3.6 




F4.5 


E4.5 


Q4.5 


F5.8 


E5.8 


Q5.8 


F8.0 




E8.0 


Q8.0 


MIPS ID 




F24 


E24 


Q24 




F70 


E70 


Q70 


12CO(1-0)'= 


13C0(2-1)'= 


Dusf^ 




1 


131.031909 


-43.736558 


263.366966 


-0.667202 


1.89903e-05 


2.54204e-06 


LI 




4.154506-06 


r\ r\r\r\r\r\ i r\r\ 

O.OOOOOc-f 00 


u 


2.11600e-05 


O.OOOOOe-l-00 


u 


2.87050e-05 




O.OOOOOe+00 


u 




0.0 


3.084e-04 


O.OOOe+00 


u 




5.220e-04 


O.OOOe-t-00 


u 


N 


N 






393 


131.096329 


-43.775259 


263.426188 


-0.654711 


1.32000C-02 


7.87000e-04 


b 




8.33412e-03 


2.59158e-04 


HI 


5.53000e-03 


3.59000e-04 


b 


3.56018e-03 




7.23856e-05 


HI 


37 


1.7 


2.543e-02 


3.492e-03 


G 




3.907e-04 


O.OOOe-fOO 


u 


N 


N 






2548 


131.183726 


-43.701613 


263.407676 


-0.559483 


1.42000e-03 


8.73000e-05 


b 




9.11296e-04 


3.42903e-05 


H2 


6.48618e-04 


5.45410e-05 


HI 


3.10650e-04 




1.79767e-05 


L2 


57 


5.1 


4.497e-02 


1.205e-03 


G 




1.165e-01 


1.191e-02 


G 


N 


N 






17590 


131.388578 


-43.868336 


263.630227 


-0.547490 


2.54057e-03 


6.21729e-05 


H2 




3.98529e-03 


1.06327e-04 


H2 


3.52000C-03 


1.78000e-04 


b 


1.29000e-03 




3.38000e-04 


b 


139 


6.0 


3.754e-01 


4.395e-02 


G 




8.675e-02 


O.OOOc+00 


u 


C 


C 


MMS2 




170299 


132.967055 


-43.592335 


264.135486 


0.511679 


5.86857e-05 


5.13993e-06 


HI 




5.54702e-05 


6.68016e-06 


HI 


2.14100e-05 


O.OOOOOe-l-00 


u 


3.13800e-05 




O.OOOOOe+OO 


u 




0.0 


2.574e-04 


O.OOOc-hOO 


u 




NaN 


O.OOOc+00 


u 


N 


N 







^Goordinatcs are given in degrees; fluxes (F) and their uncertainties (E) are given in Jy; quality flags (Q) are 
described in §[4] 

'^Angular distance in arcsec between the IRAC and MIPS centroids. 
'^Symbols are explained in §[4] when Columns 26-28 are described. 
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Table 6. VMR-D cloud IRAC catalog statistics (global) 



With good photometry 


Number of sources 


In only one IRAC band 


53461 


In only two IRAC bands 


97114 


In only three IRAC bands 


10851 


In all four IRAC bands'' 


8873 


Total number of sources 


170299 



^Including 77 recovered sources with uncertain 
photometry (see § I3.3.3p . 



Table 7. Statistics of sources toward VMR-D cloud (per band) 



Detections'^ 


3.6 /im 


4.5 /tm 


5.8 /tm 


8.0 /tm 


S/N >10 


100964 


81945 


11651 


6123 


& 24 fim 


764 


801 


634 


603 


& 24 /im & 


489 


504 


387 


356 


& 24 /xm & & Dust 


48 


50 


39 


31 


5< S/N <10 


50291 


48732 


10442 


5586 


& 24 nm 


96 


90 


68 


43 


& 24 /im & 


59 


63 


39 


22 


& 24 /im & & Dust 


3 


2 


7 


3 


Bad photometry 


19044 


39622 


148206 


158590 



Association with a 24 /tm source is within 6", with ^^CO is inside the 
5 K km s~^ co ntour level, a nd wi th a dust core is considered within the 
radii quoted in iMassi et al.l (120071 ) . 



Table 8. YSO classification of the IRAC VMR-D sources and comparison with other star forming regions. 



Region 


Area 


Distance 


YSO 


YSO 


YSO 


Class I 


fiat 


Class II 


Class III 


SFR/Area 




/I 9 \ 

{deg2) 


(pc) 


candidate 


(deg ^) 


(pc 2) 










(M0Myr pc 




U. 1 


700 
1 UU 




1 OSS 


7 'i 


an /Q 7%^ 
OU [i}. 1 /c) 


I u i 1 . o /o J 


Z4 / ^^oy.o /o j 


O/IQ /QQ 90y'\ 


0.0 


OFF ^^^f^^^r\^ 


u.uu 


700 
1 UU 


261 




9 Q 


1 1 (A 


lo ^^u.y/oj 




-LOO ^O-L.U /O ^ 


2 2 


ON-OFF'" 


0.57 


700 


372 


653 


4.4 


50 (13.3%) 


53 (14.2%) 


153 (41.1%) 


117 (31.4%) 


3.3 


ON-cloud'^ 


0.57 


700 


463 


812 


5.4 


44 (9.5%) 


61 (13.2%) 


230 (49.7%) 


128 (27.6%) 


4.1 


OFF-cloud<= 


0.60 


700 


174 


290 


1.9 


5 (2.9%) 


10 (5.7%) 


93 (53.4%) 


66 (37.9%) 


1.5 


ON-OFF'^ 


0.57 


700 


278 


522 


3.5 


39 (13.2%) 


51 (17.3%) 


142 (47.6%) 


65 (21.9%) 


2.6 


ON <^ 


0.57 


700 


465 


816 


5.5 










4.1 


OFF d 


0.60 


700 


196 


326 


2.2 










1.6 


ON-OFF ^ 


0.57 


700 


279 


490 


3.3 










2.5 


Cha II 


1.04 


178 


26 


25 


2.6 


2 (8%) 


1 ( 4%) 


19 (73%) 


4 (15%) 


0.65 


Lupus 


3.10 


150 


9 


30 


3.3 


5 (5%) 


10 (11%) 


52 (55%) 


27 (29%) 


0.83 


Perseus 


3.86 


250 


385 


100 


5.2 


87 (23%) 


42 (11%) 


225 (58%) 


31 (8%) 


1.3 


Serpens 


0.85 


260 


227 


267 


13 


36 (16%) 


23 (10%) 


140 (62%) 


28 (12%) 


3.2 


Ophiucus 


6.60 


125 


292 


44 


9.3 


35 (12%) 


47 (16%) 


176 (60%) 


34 (12%) 


2.3 



^Sources selected with the lHarvev et aL 1 1 I2OO6I . |2007|) criteria. 
''Difference counts are weighted by the corresponding solid angles. 
'^Contamination subtracted following iGutermuth et al. I l l2009l) : [Mareneo et al. I iliooi). 
"^Contamination subtracted according to lOliveira et al" 1 1 I2OO9IV 
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Table 9. VMR-D IRAC candidate YSO list: typical entries^ 



IRAC-ID 



RA(2000) 
F4.5 
E8.0 
F70 



Dec(2000) 
E4.5 
Q8.0 
E70 



Glon 
Q4.5 
MIPS ID 
Q70 



Glat 
F5.8 
M 



F3.6 
E5.8 
F24 

i2CO(l-0) 



E3.6 
Q5.8 

E24 
i3CO(2-l) 



Q3.6 
F8.0 

Q24 
Dust 



4378 131.224365 -43.926624 

9.57767E-04 3.59507E-05 



4.95515E-05 
1.199E-03 



HI 

O.OOOE+OO 



263.602173 -0.676375 1.02738E-03 2.88283E-05 HI 

HI 1.23024E-03 3.53127E-05 HI 1.36355E-03 

0.0 2.721E-04 O.OOOE-FOO u 

u -0.60 N N 



38674 131.603989 -43.704887 

3.52582E-03 6.94163E-05 



9.94560E-05 
9.515E-04 



H2 
O.OOOE-HOO 



263.599640 
H2 
261 
u 



-0.323964 2.85320E-03 7.00859E-05 



4.25391E-03 8.60992E-05 



1.1 

0.53 



4.539E-02 
C 



H2 
8.235E-04 
N 



H2 

5.84348E-03 
G 

ummsl 



'^Symbols and units as in Table [5j unless otherwise specified. Sources selected out of the working sample. 
''Spectral index of the source. 



Table 10. VMR-D MIPS candidate YSO list: typical entries. 



MIPS-ID 


RA(2000) 


Dec(2000) 


Glon 


Glat 


F24 


E24 


F70 




E70 


M 


K,-24'' 


IRAC-ID 


i2CO(l-0) 


i3CO(2-l) 


Dust 


37 


131.09596 


-43.77494 


263.42578 


-0.65472 


2.543E-02 


3.492E-03 


O.OOOE-l-00 




O.OOOE+OO 


0.0 


4.70 


393 


N 


N 




546 


131.92834 


-43.73021 


263.76633 


-0.15719 


1.163E+00 


8.047E-02 


3.519E-I-00 




3.890E-02 


1.2 


10.30 


74548 


C 





ummsll 



^Symbols and units as in Table [2l unless otherwise specified. 
^Ks-24 color index in magnitudes. 
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Table 11. Class I YSOs selected with 2MASS/IRAC/MIPS data. 





CxyZKJUKJ) 








iViir^o 






iVllr^o 


i2rin 


dust core 


IT) A 1^ 


iU 


l,aeg; 






OA T A GO 








class 


/"I n\ 
(l-Oj 




variable*^ 


17793 


131.390823 


-43.872124 


0.61 


Y 




u 


u 




C 


MMS2 


N 


18177 


131.395065 


-43.882664 


0.37 


Y 




u 


u 




C 




Y 


18652 


131.400284 


-43.891411 


0.82 


N 




u 


u 




C 




N 


19225 


131.406204 


-43.917553 


1.16 


N 




u 


u 




C 




N 


19354 


131.407394 


-43.918594 


0.66 


Y 




b 


u 




C 




N 


21173 


131.426315 


-43.456017 


0.62 


N 




u 


u 




N 




N 


21825 


131.432907 


-43.451939 


0.59 


N 


166 


G 


G 




N 




N 


22179 


131.436584 


-43.448711 


0.61 


Y 




u 


u 




N 




N 


22262 


131.137500 


-43.447865 


0.31 


Y 




u 


u 




N 




N 


23011 


131.445084 


-43.389664 


1.29 


N 


172 


G 


G 




N 




N 


38674 


131.603989 


-43.704887 


0.53 


Y 


261 


G 


u 


I 


C 


ummsl 


N 


39589 


131.613586 


-43.711151 


1.25 


N 


270 


G 


G 




C 


ummsl 


N 


41404 


131.632233 


-43.910416 


1.39 


N 




u 


u 




c 


MMS4 


N 


41429 


131.632477 


-43.549744 


1.04 


N 


290 


G 


u 




c 




N 


41786 


131.636246 


-43.911308 


0.89 


Y 




u 


u 




c 


MMS4 


N 


44510 


131.664520 


-43.380684 


0.55 


N 


316 


G 


u 




c 




Y 


47368 


131.694290 


-43.336273 


0.54 


Y 




u 


u 




c 




N 


47436 


131.694946 


-43.887310 


1.37 


N 


347 


G 


u 




c 




N 


48046 


131.701172 


-43.886124 


1.08 


N 


352 


G 


u 




c 


MMS5 


N 


50748 


131.727905 


-43.876759 


0.67 


Y 


386 


G 


u 


I 


c 




Y 


55956 


131.776840 


-43.322323 


0.50 


N 


426 


G 


u 




c 




N 


65153 


131.855423 


-43.815571 


1.03 


N 


473 


G 


G 




c 




N 


65198 


131.855774 


-43.814709 


0.73 


Y 


473 


G 


G 




c 




N 


67878 


131.877563 


-43.472893 


0.30 


N 


497 


G 


u 




c 




Y 


67896 


131.877747 


-43.467525 


0.33 


N 


498 


G 


u 




c 




N 


68154 


131.880020 


-43.898384 


1.10 


N 


503 


G 


G 




c 




N 


68231 


131.880646 


-43.896694 


0.43 


N 


503 


G 


G 




c 




N 


70412 


131.897675 


-13.356682 


0.45 


N 


521 


G 


^ 




c 




N 


72549 


131.913589 


-13.438141 


0.54 


N 




u 


u 




c 




N 


84160 


131.999054 


-43.653702 


0.43 


N 


601 


G 


u 




c 




N 


101621 


132.129349 


-43.517223 


0.85 


N 


717 


G 


u 




c 




N 


102538 


132.136902 


-43.515751 


0.79 


N 




u 


u 




c 




N 


104823 


132.155701 


-43.281319 


0.51 


N 


746 


G 


u 




c 


umms22 


N 


104849 


132.155991 


-43.280067 


0.36 


N 


746 


G 


u 




c 


umms22 


N 


106972 


132.173294 


-43.529209 


1.78 


N 


775 


G 


G 




c 


MMS9 


N 


108306 


132.184921 


-43.268906 


1.44 


N 


787 


G 


G 




c 




N 


111144 


132.209610 


-43.519905 


1.83 


N 


825 


G 


u 




c 


MMS14 


N 


112469 


132.221695 


-43.515335 


2.11 


N 


838 


G 


G 




c 


MMS16 


N 


112472 


132.221710 


-43.636738 


0.94 


Y 


840 


G 


u 


I 


c 




N 


113053 


132.227066 


-43.553127 


0.59 


Y 




u 


u 




c 




N 


114979 


132.244598 


-43.639202 


2.02 


N 


858 


G 


G 




c 


MMS17 


N 


120812 


132.302841 


-43.633789 


0.54 


Y 




b 


u 




c 




N 


120962 


132.304306 


-43.596691 


0.73 


Y 




u 


u 




c 


MMS20 


Y 


122296 


132.318268 


-43.612213 


0.35 


N 


947 


G 


u 




c 


MMS21 


N 
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Table 11 — Continued 



IRAC 


a(2000) 


5(2000) 


a 




MIPS 


Q24 


Q70 


MIPS 


12CO 


dust core'^ 


IRAC 


ID^ 


(dcg) 


(dcg) 




2MASS 


ID 






class*^ 


(1-0) 




variable'^ 


124358 


132.339371 


-44.030910 


0.73 


N 




u 


u 




c 




N 


124521 


132.340988 


-44.033684 


2.07 


N 




u 


u 




c 




Y 


137019 


132.466995 


-43.378731 


0.46 


Y 


1108 


G 


u 


I 


c 




N 


138002 


132.477051 


-43.388687 


0.87 


Y 


1120 


G 


u 


I 


c 




N 


139510 


132.491974 


-43.381367 


0.88 


N 


1130 


G 


u 




c 


umms26 


N 



^Sources selected out of the working sample. 

''Classification based on the diagram Ka vs Ks-[24] (see Fig[6ll. 

'^Dust core names are from lMagsi et aL 1 1 I2OO7I') . 

'^IRAC variability on a timescalc of six months IIGiannini et al.ll2009l'l . 



Table 12. Class I YSOs selected with 2MASS/MIPS data. 



MIPS 


a(2000) 


5(2000) 


Ks-24 


F70 


IRAC 


IRAC 


IRAC 


12CO 


dust core 


ID 


(dcg) 


(dcg) 


(mag) 




counterpart 


flags 


class'' 


(1-0) 




93 


131.29494 


-43.75684 


8.50 


N 


8642 


HI H2 HI H2 


f 


N 




130 


131.38113 


-43.89581 


8.90 


N 


16862 


b H2 b u 




N 




138 


131.38982 


-43.87410 


9.90 


N 


17713 


b H2 H2 H2 




C 


MMS2 


261 


131.60376 


-43.70510 


9.30 


N 


38674 


H2 H2 H2 H2 


I 


c 


ummsl 


308 


131.65639 


-43.87086 


10.80 


Y 


43665 


H2 b b b 




c 




315*^ 


131.66409 


-43.88627 


8.70 


Y 








c 




386 


131.72771 


-43.87690 


10.80 


N 


50748 


H2 H2 H2 H2 


I 


c 




541 


131.91966 


-43.43661 


9.20 


N 


73282 


b H2 b H2 




c 


ummslO 


546 


131.92834 


-43.73021 


10.30 


Y 


74548 


H2 H2 H2 b 




c 


ummsll 


576 


131.96498 


-43.42231 


9.50 


Y 


79572 


H2 H2 H2 b 




c 




729 


132.14119 


-43.51330 


9.50 


Y 


103109 


H2 H2 b H2 




c 




814 


132.20293 


-43.54186 


8.40 


Y 


110380 


b b b b 




c 


MMS12 


840 


132.22160 


-43.63700 


9.30 


N 


112472 


H2 H2 H2 H2 


I 


c 




884^= 


132.26060 


-44.40168 


8.80 


N 








N 




921 


132.29868 


-43.59826 


9.40 


N 


120425 


H2 H2 H2 H2 


f 


c 


MMS20 


928 


132.30367 


-43.60686 


10.80 


Y 


120963 


b b H2 H2 




c 


MMS21 


1108 


132.46693 


-43.37880 


9.30 


N 


137019 


H2 H2 H2 H2 


I 


c 




1114 


132.47073 


-43.68967 


8.90 


N 


137393 


b H2 H2 H2 









1120 


132.47690 


-43.38887 


9.30 


N 


138002 


H2 H2 H2 H2 


I 


c 




1170 


132.54562 


-43.28630 


10.00 


N 


144801 


HI H2 HI b 




c 


MMS29 



'^Classification based on the spectral index. 
''Embedded in very bright nebulosity near IRS17. 
'^Observed at 24 {rta only. 
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Table 13. Other remarkable objects selected with MIPS. 



MIPS 


a(2000) 


(5(2000) 


F24 


F70 


IRAC 


IRAC 


Distance 




dust core 


ID^ 


(deg) 


(deg) 






counterpart 


flags 


(arcsec) 


(1-0) 




31 


131, 


.03918 


-44.00486 


N 


Y 








N 




48 


131, 


,13193 


-43.83388 


Y 


Y 






3.4 


N 




146'" 


131, 


.39792 


-43.85200 


S 


Y 


18476 


b b b b 


2.8 


C 




288 


131, 


.63107 


-43.92499 


Y 


Y 


41283 


b b b b 


0.5 


C 


MMS4 


294 


131, 


.63802 


-43.92678 


Y 


Y 


41975 


b b b b 


0.9 


c 


MMS4 


300*" 


131, 


.64540 


-43.91031 


S 


Y 






6.5 


c 


MMS4 


315 


131, 


.66409 


-43.88627 


Y 


Y 






16.0 


c 




642 


132, 


.04202 


-43.35051 


Y 


Y 


90122 


H2 H2 L2 L2 


0.8 


c 




657 


132, 


.06612 


-43.78766 


Y 


Y 


93485 


b H2 b u 


2.2 


c 


ummsl6 


661=^ 


132, 


.06978 


-43.78875 


Y 


Y 


93906 


H2 H2 H2 H2 


0.9 


c 


ummsl6 


813*" 


132, 


.20084 


-42.90561 


S 


Y 






4.9 


N 




814*= 


132, 


.20293 


-43.54186 


S 


Y 


110380 


b b b b 


2.1 


c 


MMS12 


917 


132, 


.29530 


-43.49447 


Y 


Y 


120052 


H2 H2 u u 


3.2 


o 




923 


132, 


.30054 


-44.27600 


Y 


Y 






11.1 


o 




925^ 


132, 


.30151 


-43.27896 


N 


Y 


120280 


H2 H2 u u 


14.6 


c 




951 


132, 


.32057 


-43.93255 


Y 


Y 


122524 


u b u L2 


0.1 


c 




991*" 


132, 


.35916 


-43.28617 


S 


Y 


126339 


b b b b 


1.5 


c 


MMS22 


998 


132, 


.36937 


-44.07478 


Y 


Y 


127325 


b b b b 


1.5 


c 


MMS23 


1018*" 


132, 


.38667 


-44.18056 


S 


Y 






3.6 


c 


MMS26 


1247° 


132, 


.65999 


-43.82816 


N 


Y 


156273 


u b H2 H2 


12.0 


N 




1306 


132, 


.87169 


-43.52266 


Y 


Y 






9.3 


N 





^AU the sources fall in the area mapped with IRAC and MIPS (24/jm and 70/im). 

'"So urces saturated at 24 fim: # 146, 300, 813, 814, 991, 1018 corresponding to IRS 16, 17, 18, 19, 20, 21 of lLiseau et~all 
l|l992h . respectively. 

'^Classified as Class II with IRAC. 

"^Barely visible at 70 /im, in the other bands this position corresponds to the boundary of a spherical blob ( 2' in 
diameter); the IRAC source could not be a genuine association due to the large distance (column 8). 

"Barely visible at 70 ^m; at 24 fim there is a very strong diffuse emission and a faint object is visible at a distance 
of ~ 12" NW. This faint source is excluded from the 24 /tm list because of the constraints imposed to the PSF fitting 
parameters (see ij 13.51 1. 



